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Recently a cluster Monte Carlo algorithm has been used very successfully in the two- 
dimensional Edwards-Anderson (EA) model. We show that this algorithm and a variant 
thereof can also be used successfully in models with a non-zero spin glass transition temper- 
ature. The application of such algorithms to the site-diluted EA model in three dimensions is 
discussed and the efficiency of the two algorithms is compared among each other and to par- 
allel tempering. Finally, we give evidence for a spin glass transition in the three-dimensional 
site-diluted EA model with Gaussian couplings at a site occupation of p — 62.5%. 

§1. Introduction 

The determination of equilibrium properties of finite dimensional spin glass mod- 
els relies heavily on Monte Carlo (MC) simulations. However, these simulations have 
shown to be a very difficult task due to the strong critical slowing down at the phase 
transition and the existence of big free energy barriers in the low temperature phase. 
There are few algorithms that sample the relevant regions of the phase space in 
spin glass models efficiently enough to allow for the simulation of systems with sizes 
around 16 3 to 20 3 spins with reasonable statistics. 

The use of cluster algorithms which revolutionized the simulations in ordinary 
spin models^^ hasn't shown to be very effective in spin glass models due to frus- 
tration and the lack of simple symmetries. Recently, however, Houdayer intro- 
duced a cluster algorithm combined with the exchange MC method (also known 
as parallel tempering (PT)) which has shown to be extremely efficient in the two- 
dimensional Edwards- Anderson (EA) modelf^*) The main idea of this algorithm is 
already present in 03) ,|f)J) and is to simulate two or more replicas of the same system 
and to use their mutual overlap to generate clusters. The main drawback mentioned 
in 3) is that the algorithm performs poorly for the three-dimensional EA model be- 
cause of a site percolation problem and therefore the method seemed to work only 
for systems with a spin glass transition at zero temperature. 

In this paper we show that Houdayer 's algorithm and a related algorithm de- 
scribed below do also work efficiently in a certain class of spin glass models with a 
finite spin glass transition temperature. Here we will focus on the three-dimensional 
site-diluted EA model given by the HamiltoniaiP 

H = - ^ .J, j( ,< :rr,rr : r (1-1) 

<ij> 



*' A slightly modified version of the original Swendsen-Wang replica cluster algorithm has even 
shorter autocorrelation times than Houdayer's algorithm in the two-dimensional EA modelP' 
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where the sum is over all the nearest neighbor sites of a simple cubic three-dimensional 
lattice of length L, the o% = ±1 are Ising spins and the ej are the site occupation 
variables, i.e. ej = for an empty and ej = 1 for an occupied site. For the couplings 
Jij we use a Gaussian distribution with mean zero and standard deviation unity. The 
fraction of occupied sites is denoted by p. At this point we would like to mention 
that there are other diluted models for which such replica cluster (RC) algorithms 
work very efficiently as e.g. the Viana-Bray mode^ 8 ' with fixed connectivity z = 3 
and also the link-diluted EA model in three dimensions PES 

The site-diluted EA model in three dimensions is clearly a model of physical 
relevance*^ . One prototype of a real spin glass consists of a small fraction of randomly 
diluted magnetic moments dispersed in a non-magnetic metal and hence site disorder 
is an inherent property of such a system. It is therefore desirable to check numerically 
whether the EA model is stable with respect to site-dilution, as one might expect 
on general groundsPJ Note that this is in contrast to the expectations in the site- 
diluted Ising model, where site-dilution shows to be a relevant perturbation that 
changes the critical behavior of the system.^ In the following we are going to 
discuss a variation of Houdayer's RC algorithm used in this study. We show that 
it is very efficient in thermalizing the spin glass model defined through Eq. 
compare it to the original version of Houdayer's algorithm^ and to PtP^J Finally, 
we discuss possibilities and limits of this kind of RC algorithms and show results for 
the finite-size scaling of the correlation length. 



§2. Algorithm and its limitations 



The idea of the algorithm is as follows. Two independent replicas of the system 
are simulated simultaneously. Hence each site of the lattice is associated with two 
spins, one from each replica, and can be in one of four spin states (++), (+— ), 

( — h), and ( ). We call the sites where the replicas disagree, (H — ) and ( — h), 

active sites. In the case of a non-zero external magnetic field clusters of active sites 
are constructed and flipped on both replicas. Such cluster flips interchange (H — ) 
with ( — h) and vice versa. In the case of no external magnetic field also the sites 
where the replicas agree (non-active sites) can be used to construct clusters and we 
will comment on this possibility below. In order to ensure ergodicity of the algorithm 
we update each replica independently by wrapping the cluster moves within a PT 
as proposed by HoudayerP 

The clusters themselves are formed as described in|HJ), i.e. they are Wolff- type 
clusters grown on the active sites of the lattice. However, having a disordered Hamil- 
tonian we need to generalize the procedure slightly as follows. The clusters are grown 
on the active sites by adding links with probability p a dd- 

fl-exp(4/3£y) if Eij<0 
Padd = \0 else, (2 " 1} 

where EL- = —Jij€iCjaiaj is the energy contribution of the given link and (5 is the 
*' The EA model in the original worlPl is actually formulated as a site disordered model. 
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inverse temperature. This is the only difference to Houdayer's algorithm, where the 
links between two active sites are always added to the cluster. The clusters we have 
defined through Eq. (|2-lj) flip freely i.e. the cluster move is always accepted. Note 
that also more general clusters might be constructed along the lines of I14|) . where 
the size of the clusters can be controlled in a very general manner. In the absence 
of an external field one may as well grow clusters on the non-active sites, as the 
constraint to the active sites is just to ensure that the total magnetic energy of the 
two replicas is conserved, as (H — ) goes to ( — h) and vice versa. Without external 
field one may hence also do a full Swendsen-Wang updating of all clusters, which in 
diluted systems may be helpful to update clusters of all sizes 

The algorithm we proposed is a valid MC procedure for general Ising spin sys- 
tems. However, for frustrated interactions it will generally not perform well when 
the site percolation threshold psp of the underlying interaction graph is below 0.5, 
because in this case only very small (or equivalently very large) clusters are formed. 
If, on the other hand, psp > 0.5 clusters of all sizes can be formed in the disordered 
phase, allowing for a fast decorrelation of the configurations. 

§3. Results and Conclusions 

For the three-dimensional site-diluted EA model the RC algorithms are perform- 
ing clearly better in thermalizing a configuration than normal PT, as illustrated in 
the left panel of Fig. ^ f° r a 12 3 system with p = 0.5 at roughly 0.6T C . Although 
one MC step in a RC algorithm is more expensive than for PT, the performance in 
absolute computer time is by far better for the RC algorithm. The right panel of 
Fig. ^ shows how the cluster algorithm defined through Eq. (|2-H) compares to Hou- 
dayer's algorithm. The equilibration of the energy on a 16 3 system with p = 0.625 
at roughly 0.8T C is significantly faster for the Wolff-type clusters. 



<M rwx> 




Fig. 1. Equilibration of the energy. Left panel: Comparison of the RC algorithm with Wolff- type 
clusters (open squares) with PT (filled circles). Right panel: Comparison of the RC algorithm 
with Wolff-type clusters (open squares) with Houdayer's clusters (filled squares). Note that the 
two comparisons are done on different system sizes. 

Finally, Fig. |2] shows the finite size behavior of the overlap correlation length 
£l/L, as defined in ll5f) . for system sizes ranging from 4 3 to 20 3 at a site dilution of 
p = 0.625. For small systems the curves do not cross at a unique point due to finite 
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size effects. For larger systems, however, the curves start to cross very close to each 
other, such that this data presents clear evidence for a phase transition at T c ~ 0.51. 
The value of at T c is in good agreement with the simulations of the undiluted 
three-dimensional ±«^^ and Gaussiarf^ EA model. More details are given inll8|). 

RC algorithms work efficiently in some 
systems with a non-zero spin glass temper- 
ature, typically in models with a low con- 
nectivity. This opens a new perspective 
to investigate spin glass models that have 
not been well accessible by more traditional 
MC methods, like site- or bond-diluted EA 
models in three dimensions. Moreover, we 
have shown that using a Wolff- type defini- 
tion of the clusters can have a benefit with 
respect to Houdayer's definition. A more 
precise determination of the critical prop- 
erties of the three-dimensional site-diluted 
EA model will allow to give clearer results 
about its universality class. Finally, there 
remains the question whether the fact that 
efficient cluster algorithms for some very 
specific spin glass models exist is just an incident or not. 
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Fig. 2. Data for the correlation length 
£ L /L for p = 0.625. 
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